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ABSTRACT 

We present 21-cm observations and models of the neutral hydrogen in NGC 4565, a nearby, edge-on 
spiral galaxy, as part of the Westerbork Hydrogen Accretion in LOcal GAlaxieS (HALOGAS) survey. 
These models provide insight concerning both the morphology and kinematics of H I above, as well 
as within, the disk. NGC 4565 exhibits a distinctly warped and asymmetric disk with a flaring layer. 
Our modeling provides no evidence for a massive, extended H I halo. We see evidence for a bar 
and associated radial motions. Additionally, there are indications of radial motions within the disk, 
possibly associated with a ring of higher density. We see a substantial decrease in rotational velocity 
with height above the plane of the disk (a lag) of — 40 1 20 km s _1 kpc -1 and — 30 km s -1 kpc -1 
in the approaching and receding halves, respectively. This lag is only seen within the inner ~4.75' (14.9 
kpc) on the approaching half and ~4.25' (13.4 kpc) on the receding, making this a radially shallowing 
lag, which is now seen in the H I layers of several galaxies. When comparing results for NGC 4565 
and those for other galaxies, there are tentative indications of high star formation rate per unit area 
being associated with the presence of a halo. Finally, H I is found in two companion galaxies, one of 
which is clearly interacting with NGC 4565. 

Subject headings: galaxies: spiral - galaxies: halos - galaxies: kinematics and dynamics - galaxies: 
structure - galaxies: individual (NGC 4565) - galaxies: ISM 



1. INTRODUCTION 

Understanding the origins and evolution of extra- 
planar gas is key to interpreting how spiral galaxies are 
affected by their environments. Whether galaxies ex- 
ist in dense clusters or groups, or are relatively isolated, 
extra-planar gas is the bridge between galaxy disks and 
the intergalactic medium (IGM). Is this gas ejected from 
the plane of the disk, driven by star formation as de- 
scribed in the galactic fountai n dShapiro fc Fieldl I1976L 
lBregmar] ri980) and chimney (jNorman fc Ikeuchi I 1989ft 
models? Does some of it result from accret ion from 
external sources (e.g. iKeres fc Hernquisll 120091 ) such as 
companion galaxies or the IGM itself? If so, how would 
these internal and external components interact with 
each other? To answer these questions, a thorough study 
of extra-planar gas morphology and kinematics in a large 
number of galaxies is necessary. 

Extra-planar c ompo nents including hot gas (e.g. 
Tull mann est al.1 12006ft . r elativistic particles (e.g. 
Irwin et al.lll999ft . dust Ce g. IHowk fc Savagdll999ft. ion- 
ized hydrogen (e.g. lRandHl996l; iRossa fc Dettmarl 12003ft 
as well as neutral hydrog en (H il (e.g. iSwaters et al.l 
119971 lOosterloo et al.N2007l ) have been observed. For all 
except H I there is a correlation between their presence 
and star formation in the disk, both in localized regions 
as well as globally. This indicates disk-halo flows akin 
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to those in the aforementioned galactic fountain and 
chimney models as likely origins. However, such trends 
have not been reported for Hi, and in fact there is 
evidence that at least some ext ra-planar H I has an ex- 
terna l origin in the Milky Way (IWakker fc van Woerdenl 
11997ft and nearby galaxies (jSancisi et al.l 12008ft . This is 
an issue that motivates this work. 

In addition to probing the origins of different halo com- 
ponents through their morphology, we may also gain in- 
formation through understanding their kinematics. In 
galactic fountain-type models, gas is launched from the 
plane of the disk and moves to larger radii. In order to 
conserve angular momentum, there must be a decrease in 
rotational velocity with height. This vertical gradient in 
rotational velocity is referred to as a lag. However, other 
physical processes may affect the magnitude of such a 
lag (see below); hence, measurement of lags should help 
to constrain the physics of extra-planar components, and 
may shed light on important elements of galaxy evolu- 
tion. 

IHeald et all (|2007ft measured lags in three extra-planar 
Diffuse Ionized Gas (DIG) layers, resulting in a trend 
of low star formation rates (SFR) with steeper lags. If 
H I and DI G laye rs share similar kinematics and the 
IHeald et al.1 (|2007ft trend holds, then, the same should 
be expected for H I. 

Simple dynamical models of disk- halo cycling, where 
only ballistic effects are c onsidered (|Collins et al.l 120021 
iFraternali fc Bmrievl 12006ft . have produced lags signifi- 
cantly shallower than those observed. These discrepan- 
cies between models and observations indicate that addi- 
tional factors must be considered. Firstly, internal com- 
plications such as those produced by extra-planar pres- 
sure gradients or magnetic tension (|Benjaminl 12002) have 
yet to be examined thoroughly, but could potentially ex- 
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plain these discrepancies without resorting to external 
factors. 

Explanations involving infalling gas such as interaction 
of disk-halo cycled gas with a possible hot, low angular 
momentum, low metallicity halo predi cted by cosmologi- 
cal simulations (jMarinacci et al.ll2011h and the hydrody- 
namic simulation of disk formation in iKaufmann et all 
(2006) have successfully reproduced the observed lag in 
NGC 891. 

However, while much attention has been focused on re- 
sults for NGC 891, o bserved lags differ a mong the three 
galaxies considered in iHeald et al.l (|2007| ) and one would 
like to understand the physical cause of such variations. 
More relevant here, extensive modeling of extra-planar 
H I kinematics has until recently been carried out for 
only a small number of edge-on galaxies, using observa- 
tions with a rang e of resolutions and sen sitivities. These 
include NGC 891 (jOosterloo et alJ2Ml . where DIG and 
H I lags agree, pointing to similar origins for both compo- 
nents, and other galaxies summarized in § 16.41 Further- 
mo re, the trend of st eeper lags with lower SFRs found 
bv IHeald et al.1 (|2007f ) needs to be further tested, while 
variation of lags within extra-planar gas layers, such as 
radial gradients, may also constrain their origin. With 
the aforementioned evidence for an external source for 
at least some extra-planar H I, it is imperative to under- 
stand how kinematic information may constrain these 
two possible origins of extra-planar gas. 

The Westerbork Hydro gen Accretion in L Ocal GAlax- 
ies (HALOGAS) survey ijHeald et all 120111 ). which tar- 
gets 22 edge-on or moderately inclined spiral galaxies for 
deep, uniform 10x12 hour observations using the West- 
erbork Synthesis Radio Telescope (WSRT), will greatly 
increase this number. This survey will also allow us to 
establish whether there is any connection between extra- 
planar H I and star formation, as well as the degree to 
which contributions from external origins are relevant. A 
primary goal of HALOGAS is to estimate the rate of H I 
accretion in spiral galaxies. 

Aside from the characterization of extra-planar gas, 
one goal of HALOGAS is to investigate morphological 
features such as bars and rings as well as their kinematics. 
These features will be discussed throughout this paper. 

1.1. NGC 4565 

NGC 4565 is cl assifie d as a SAb galaxy 
(|de Vaucouleurs et al.l pED) with a SF R of 0.67 
M Q yr" 1 calculated in IHeald et al.l (|2012| ) using total 
infrared (TIR) flux measurements. This places it near 
the low star-forming end of the HALOGAS sample. 
However, based on its rotation curve, it is among the 
most massive, which would presumably hinder gas being 
ejected to large heights above the disk. This combination 
gives it substantial importance in determining overall 
trends involving extra-planar H I and star format ion. 

At a distance of 10.8 Mpc (|Heald et al.ll20rl . NGC 
4565 is relatively nearby, and among the best stud- 
ied edge-on spirals. It has been observed in opti- 
cal continuum (e.g . iKormendv fe Barentind I2010T). Ha 
(I Rand et al.1 11992b . rad io continuum ([Kodilkar et al.l 
12008]) CO line emission (jNeininger et al.lll996D . and X- 
rays (jVogler et al .111 9961 ). 

In the optical, NGC 45 65 is seen to have a box-shape d 
bulge indicating a bar (jKormendv fe Barentind 120101) . 



Additionally, iNeininger et al.l ([1996) found radial mo- 
tions in CO emission along the minor axis, which they 
suggested could be due to a bar or a spiral density wave. 
We see evidence for the bar in H I as will be discussed in 

§E|__ 

iKodilkar et al.l (|2008| ) detect a ra dio continuum halo 
extending up to 3 kpc. In contrast, I Rand et all (|1992f ) 
found that NGC 4565 was lacking a smooth DIG halo, 
which is consistent with its low star formation rate. 

Unexpectedly, substantial X-ray emission, generally 
not seen in galaxies without enhanced star formation 
activity, was detec ted both within and above the disk 
(jVogler et al. 1996). The nature of this X -ray emission 
is still unknown. 

Finally, H I observations of NGC 4565 were previ- 
ously performed using the Very Large Array (VLA), 
but detailed models were not created from those data 
(jRupenl ri991h . However, a comprehensive description 
was provided in that paper, mentioning several fea- 
tures we now see in greater detail, and with a higher 
level of confidence. Additional WSRT data were taken 
(iDahlem et al.l 12005ft. w hich were further interpreted in 
Ivan der Hulst fc Sancisil (|2005l ). These WSRT data are 
supplemented with the HALOGAS data to produce the 
final dataset presented here. Q] summarizes parameters 
for NGC 4565. 

2. OBSERVATIONS AND DATA REDUCTION 

Here we provide a brief explanation of the observa- 
tions and data reduction process. A detailed description, 
which a pplies to all of the HALOGAS galaxies, may be 
found in IHeald et all (|201lD . 

NGC 4565 was observed using the Westerbork Syn - 
thesis Radio Telescope (WSRT) bv IDahlem et al.l (2005) 
as well as through the HALOGAS survey. The archival 
observations were done in the Traditional WSRT con- 
figuration, using four distinct spacings between the fixed 
and movable antennas, starting at 36 m and incremented 
by 18 m. Observations obtained as part of HALOGAS 
were in the Maxi-short configuration with baselines rang- 
ing from 36 m to 2.7 km in order to maximize sensitivity 
to faint, extended emission. Of the 10 fixed antennas, 

9 were used and spaced at 144 m intervals on a regular 
grid. The total bandwidth is 10 MHz with 1024 channels 
and two linear polarizations. The total observing time is 

10 x 12 hours, and observing dates are listed in Table El 
Da ta reduction was performed in Miriad (|Sault et al.l 

1995). Ima ges were created with a variety of weighting 
scheme s, creating m ultiple data cubes. Clark deconvo- 
lution (|ClaIk1 [l980h was performed using mask regions 
defined on the basis of unsmoothed versions of the image 
cube. Offline Hanning smoothing led to a final velocity 
resolution of 4.12 km s _1 in a single channel. The la 
rms noise in a single channel of the full resolution cube 
is 0.31 mjy bm _1 , corresponding to a column density of 
N H i = 5.1 x 10 18 cm" 2 . 

To optimize our resolution, we use a cube with a 
robust parameter of —2. The beam for this cube is 
24.27x11.45"; 1" = 52.4 pc (1270x600 pc) with a po- 
sition angle of —0.47°. By using this cube, we sacrifice 
some sensitivity to faint extended emission. However, 
we later examine and apply the same models to a second 
cube where a Gaussian uv taper yielding 44 x 34" resolu- 
tion is used. Primary beam correction has been applied 
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TABLE 1 

NGC 4565 Parameters 



Parameter 



Value 



Reference 



Distance (Mpc) 10.8 a 

Systemic velocity (km s — 1 ) 1220 

Inclination 87.5 ° 

SFR (M Q yr- 1 ) 0.67 

Morphological Type SAb b 

Kinematic Center a (J2000.0) 12h 36m 21s 

Kinematic Center 8 (J2000.0) 25d 59m 10s 

D 25 (kpc) 16.2 

Total Atomic Gas Mass (1O 9 M ) 9.9 c 



Heald et al. (2011) 
This work 
This work 
Heald et al. (2011) 
dc Vaucouleurs et al. (1991) 
This work 
This work 
dc Vaucouleurs et al. (1991) 
This work 



a Distance is the median value of distances found on the NED database, excluding those 
obtained using the Tully-Fishcr relation. 

k The A classification indicating no bar is likely erroneous based on Neiningcr et al. 
1996, Kormcndy &z Barcntinc 2010, and this work. 
c Includes neutral He via a multiplying factor of 1.36. 



TABLE 2 

Observational and Instrumental Parameters 



Parameter 



Value 



Observation Dates - HALOGAS 



Observation Dates — Dahlcm 



Pointing Center 

Number of channels 
Velocity Resolution 
Full Resolution Beam Size 

RMS Noise - 1 Channel (24.27x11.5") 
RMS Noise - 1 Channel (44x34") 



2009 Mar 10-11 
2009 Mar 11-12 
2009 Mar 15-16 
2009 Apr 1-2 
2009 Apr 3-4 
2009 Apr 12-13 
2009 Apr 24 
2009 Apr 27-28 
2003 Feb 16-17 
2003 Apr 14-15 
2003 May 14-15 
2003 May 19-20 
12h 34m 21.071s 
25d 59m 13.50s 
148 

4.12 km s" 1 
24.27x11.5" 
1270x600 pc 
0.31 mjy bm 
0.24 mjy bm -1 



to both cubes using the Miriad task linmos. In the case 
of mome nt maps, wh i ch ar e created using methods de- 
scribed in lSerra et a l. ( 2012), this is done in the last step. 

3. THE DATA 

Upon initial inspection of the zeroth moment map 
(Figure [1]), H I emission can consistently (i.e. along most 
of the disk) be seen up to heights of 1.5-1.75' in the full- 
resolution and smoothed cubes respectively. As will be 
shown in § dj this does not mean that an extended H I 
halo is necessary to fit the data. Also evident in the 
zeroth moment map is an asymmetric warp component 
across the line of sight. This warp is most prevalent in 
the receding (NW) half of the galaxy, where the disk 
tilts upward and then flattens at large radii. In the faint 
emission, a connection is seen between NGC 4565 and the 
companion in the northeast quadrant (IC 3571), which 
was also noted by Ivan der Hulst fc Sancisil ((2005). This 
interaction may contribute to the asymmetry of the warp. 

The warp across the line of sight is also clear in 
the channel maps shown in Figure [2j Also seen more 
clearly in the channel maps is a substantial asymmetry 




50" 30 s 36 m lCT 13"35 m 50 s 35 m 3CT 

R.A. (2000.0) 

Fig. 1. — The Hi zeroth-moment map displaying the full- 
resolution (black contours and grayscale) and smoothed 44" X 34" 
(gray contour) cubes. Full resolution contours begin at 3.0 X 10 19 
cm -2 and increase by factors of 2. The smoothed contour is at a 
level of 1.1 X 10 19 cm -2 . Both beams are shown in the lower left- 
hand corner. Two companion galaxies may be seen, NGC 4562 in 
the southwest corner and IC 3571 north of the center. Note the 
strong, asymmetric warp across the line of sight as well as the con- 
nection between NGC 4565 and IC 3571. Diagonal lines parallel 
to the minor axis indicate slice locations shown in Figure [3) 

in the distribution of gas in the approaching and receding 
halves. 

There exist two nearby companions within the field of 
view: IC 3571 near the northeast quadrant of the galaxy 
with a velocity range of 1225-1299 km s _1 and NGC 4562 
in the southwest with a velocity range of 1266-1427 km 

There is a large velocity spread seen in minor axis 
position-velocity (bv) diagrams near the center of the 
galaxy [Figure [3] (column 1, row 5)]. As we will show in 
§ I4.1.3) this velocity spread is kinematic evidence for a 
bar with an orientation that is at least partially along 
the line of sight with radial motions. 

There also exist indications of radial motions not asso- 
ciated with the bar. These are best seen as a slanting of 
the data contours on the systemic side in most panels of 
Figure [3l The direction of the slant at minor axis offsets 
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Fig. 2. — Representative channel maps of the data with a res- 
olution of 24.3" Xll. 5". Contours begin at 2<r (0.61 mjy bm ) 
and increase by factors of 2. The —2a contour is also shown with 
dashed lines. The velocities in km s — 1 are given in each panel and 
the velocity spacing is approximately 7 resolution elements. Note 
the slanted outer edges present in both halves, which correspond 
to the projection of the warp perpendicular to the line of sight as 
well as the asymmetric gas distribution in each half. The beam is 
shown in the bottom left panel. 

< 40" is different from that at minor axis offsets 40" - 
80" . Our modeling will show that this indicates a sign 
change in the direction of radial motions at larger radii 
rather than with height above the plane. 

Crucial to certain aspects of the modeling, such as ra- 
dial motions, is the distinction between the near and far 
side of NGC 4565. Fortunately, from the orientation of 
the dust lane seen in the 2MASS Large Galaxy Atlas 
(|Jarrett et al.ll20 03) , we can infer that the side northeast 
from the plane of the disk is nearer than the southwest- 
ern side. This side corresponds to the positive side along 
the minor axis in the figures. 

Modeling of the features listed above will be discussed 
further in the next section. 

Finally, from our continuum data we simply present 
a vertical profile in Figure [4 The detectab l e ext ent is 
comparable to that found bv lKodilkar et al.l (2008). 



4. THE MODELS 



Models were c reated using the t ilted ring fitting soft- 
ware (TiRiFiC) (|J6zsa et alj|200 7H which allows for a 

6 http:/ /www. astron.nl/~jozsa/tirific/ 



X 2 minimization fit. Aside from its fitting capabilities, 
this software employs methods simil ar to the Groningen 
Imag e Processing System (GIPSY) ([van der Hulst et all 
1992) task galmod. For each ring, parameters such as 
the central position, systemic velocity, inclination, po- 
sition angle, surface brightness, rotational velocity, ra- 
dial velocity, velocity dispersion and scale height may be 
specified. 

Contrary to models created using galmod, those cre- 
ated using TiRiFiC need not be axisymmetric. TiRiFiC 
allows the user to divide the model galaxy into wedges 
of azimuth in order to model asymmetries throughout 
the disk without requiring the creation of multiple mod- 
els. Additionally, TiRiFiC allows for a gradient in the 
rotational velocity with height. TiRiFiC can also super- 
impose gaussian components onto the disk, which can be 
used to resemble bars and spiral structure in the model. 
Naturally, as is the case for tilted ring models in general, 
these are greatly oversimplified and should only be used 
as tools to approximate the morphology of these features 
and not to fully explain their kinematics. 

Although a majority of the modeling of NGC 4565 
was done in TiRiFiC, the earliest modeling steps were 
done using standard GIPSY tasks. This is only because 
TiRiFiC was still undergoing substantial improvements 
at the beginning of the modeling process. 

An initial estimat e for the systemic velocity was taken 
from iRupenl (|1991h . The rotation curves were deter- 
mined via inspection of the terminal side of major axis 
position- vel ocity (lv) diagrams, akin to the envelope trac - 
ing method (|Sancisi fc AllerjH979L ISofue fc Rubinl[200ll) . 
Surface brightness estimates were determined using the 
GIPSY task radial. The run of the position angle was 
initially estimated by eye and inclinations close to 90 ° 
were considered. Initial estimates for the exponential 
scale height were determined solely by a fit to the vertical 
profile (Figure [5]) . All of these parameters were later re- 
fined both by comparing observed and modeled position- 
velocity diagrams, channel maps, zeroth-moment maps 
and vertical profiles. 

TiRiFiC ' s automated fitting capabilities were used 
when appropriate. However, at the level of detailed anal- 
ysis of lopsided and asymmetric disks, which may be 
described by including non-standard parameterizations, 
this possibility is useful only to a certain degree. Nev- 
ertheless, if provided with reasonable initial constraints, 
TiRiFiC may expedite the fitting process. The degree to 
which additional refinement by eye is necessary depends 
greatly on any distinct features present in the region of 
fitting. 

4.1. Individual Models 

Several basic models are initially considered. These in- 
clude a simple one-component (i.e. single disk with a con- 
stant exponential scale height and no warp or flare), flare, 
a substantial warp component along the line of sight, and 
a model with a second, thickened disk simulating an ex- 
tended global halo. The one-component and flare models 
are labeled as 1C and F in figures. 

Upon examination of the zeroth moment map (Fig- 
ure [TJ) as well as the vertical profile (Figure [S]), the latter 
two models (not shown) may be eliminated quickly due 
to a lack of thickening of the gas layer at small to inter- 
mediate radii that is a feature of these models as well as 
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Fig. 3. — Position-velocity diagrams parallel to the minor axis showing the data, one-component (1C) model, and various flaring (F) 
models. The Fi model is the flare model with disk inflow at small to moderate radii beyond the bar, and the Fio model is the flare model 
with both inflow in the inner parts of the disk and outflow in at large radii. The second to last column shows the Fio model with a global 
lag of —10 km s _1 kpc — 1 , while t he fin al column shows the Fio model with an optimize d radia lly varying lag. Details for models involving 
radial motions are described in § 14.1.51 and details of the lag models are described in § 14.1.61 All models shown here include a bar with 
an orientation 30° from the line of sight and a radial motion of —20 km s —1 . Contours are as in Figure [2] Slice locations are given in 
the first column and can be seen in Figure [T] Arrows indicate the direction emission is moved due to the addition of radial motions. The 
approaching half is SE while the positive minor axis offset is NE as seen in Figure [T] 
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Fig. 4. — The vertical profile of the Hi continuum at 30" resolu- 
tion corresponding to major axis values between —4' and +4'. The 
FWHM is approximately 1.2' or 3.8 kpc. 
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Fig. 5. — The vertical profile of the data, one-component (1C), 
and flaring (F) models corresponding to major axis values between 
—4' and +4'. Especially note the fit of each model near the peak 
of the profile as well as near the wings. 

poor fits to position-velocity diagrams and channel maps. 
The observed vertical profile is almost entirely devoid of 
wing structure, making it difficult to fit while maintain- 
ing the correct curvature and width near the center with 
such models. A substantial warp component along the 
line of sight places flux high above the midplane on the 
systemic side of the bv diagrams as is seen in the data. 
However, it also creates indentations along the major axis 
at large radii in channel maps as well as on the systemic 
side of bv diagrams. These are not seen in the data. As 
for a halo model, adding a second thicker component not 
only increases the width of bv diagrams on the systemic 
side, but also increases the width closer to the terminal 
side. Again, this is not seen in the data. For these rea- 
sons, only the 1C and F models are considered further 
and described below. 
Although multiple types of models are considered, 



some parameters remain fixed. These include the cen- 
tral position (12h 36m 20.63s, 25d 59m 16.8s), systemic 
velocity (1220 km s _1 ), velocity dispersion (10 km s _1 ), 
the run of the position angle and inclination, the surface 
brightness profile, and the rotation curve. All models in- 
clude a bar and a small-scale feature that are discussed 
in 5 HXB1 and jXj 

Due to asymmetries, the disk is divided into two halves 
along the minor axis, and each half independently mod- 
eled in what follows. 

4.1.1. One Component Model 

A model consisting of a single component with an ex- 
ponential scale height of 470 pc was created. This model 
is a reasonable approximation to the data, but with de- 
ficiencies. The first of these may be seen in the vertical 
profile (Figure [SJ where the 1C model is too flat near the 
top compared to the data. If the scale height is decreased 
to match the data near the center, then the model be- 
comes too narrow in the wings. No other combination of 
scale height and inclination improves the fit. In the bv 
diagrams shown in Figure 02 the 1C model is too thick 
at terminal and intermediate velocities when compared 
to the data. Furthermore, an excess of emission in the 
1C model can be seen at major axis offsets within ±5' 
of the center offsets of the lv diagrams corresponding to 
-16" and -32" in Figure l6al and b. Finally, in Figure[7j 
one may note that, while the minor axis thickness is well 
matched at velocities close to systemic, the 1C model is 
too thick at small and intermediate major axis offsets. As 
will be shown, these issues are all remedied by allowing 
the scale height to vary in the form of a flare. 

4.1.2. Flaring Model 

Instead of using a single exponential scale height 
throughout the disk, if we allow the disk to flare [scale 
height = 4" ( 200 pc) near the center, increasing to scale 
height = 12" ( 600 pc) in outer radii], the issues men- 
tioned above are all but eliminated, as may be seen in 
the third columns of Figures IB"al and b. Other variations 
in the scale height were tested, but this was found to be 
the best fit to the data. Figure 03 shows the radial varia- 
tion of scale height used in our flaring models. Note that 
the scale height increases in steps rather than a straight 
line. A smoothed version would produce nearly iden- 
tical results, but arbitrarily smoothing the scale height 
distribution could imply that the models are better con- 
strained than they are in reality. Thus, we include the 
steps as they are a better representation of the modeling 
process and limitations. Nevertheless, the scale height is 
seen to rise in a nearly linear fashion. 

Although some improvements with the flare model are 
seen in lv diagrams where the excess emission in the 1C 
model is not seen, the most convincing evidence for im- 
provement is also seen in bv diagrams and channel maps 
(Figures [3] and [7]). In Figure 03 one may note a widening 
on the systemic relative to the terminal side of the data 
best seen in panels within ±3', which is significantly bet- 
ter matched by adding a flare. Additionally, in Figure UJ 
it is seen that there is a thickening at large major axis 
offsets in the channel maps that is best reproduced via 
the addition of the flare (e.g. the panel corresponding to 
1381 km s- 1 ). 
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Fig. 6a. — Position-velocity diagrams of the approaching (A) and receding halves (B) parallel to the major axis showing the data, one- 
component (1C), and flaring (F) mo dels. The final t wo columns (Fio and Fio + Lag) display the flare model with radial motions and then 
with the optimal lag described in § 14.1.51 and 14.1.61 All models shown here include a bar with an orientation 30 ° from the line of sight 
and a radial motion of —20 km s . Contours are as in Figure [2] The heights above the plane of the disk are given in the first column. 
Note that the lag i s di fficult to detect in these diagrams, primarily due to its being only within the inner 3-4'. More convincing evidence 
is shown in Figure 1131 

Key parameters for the models described above, with 
the 1C and F models differing only in scale height be- 
havior, are shown in Figure [5] 

Up to this point, we have discussed modeling only 
global features found throughout the disk. Now we will 
discuss the modeling of more localized features such as 
a bar and adding further refinements such as radial mo- 
tions. Given that Figures [3] through 5 provide the most 
detailed representation of the models, these figures all 
include the features that will be described below. Ab- 
breviated representations of the models excluding them 
will be shown in order to show their effects. 

4.1.3. Including a Bar 

As mentioned in § [T] and [3J there are indications of 
a bar in NGC 4565, which we will show can also ex- 
plain the large velocity spread in panels near the cen- 
ter in Figures [3] and [TO] The effects of adding a bar 
are most clearly seen in the latter Figure, which shows 
no-bar and bar models side-by-side. We model this bar 
by adding a series of ID gaussian emission components 
starting at the center and moving outward radially along 



some azimuthal direction in the disk. The dispersion of 
the gaussian component is specified along the azimuthal 
angle and is 15" in our models. The amplitude is also 
specified. Additionally, parameters which are specified in 
rings as described above may also be set, notably the az- 
imuthal position and radial velocity. It should be noted 
that these parameters allow for an approximation of a 
bar, and do not account for more complex streaming 
motions or morphology. Figure [5] shows how this bar 
appears in a face-on view of our best model. 

We examine several orientations with re- 
spect to the line of sight, keeping in mind that 
iKormendv k. Barentind (|2010l ) assume a line of sight 
(i.e. end-on) orie ntation due to the bo xy bulge seen in 
the optical, while lAthanassoulal (|2005f) states that such 
boxy bulges are indicative of components across the 
line of sight (i.e. side-on), but possibly only as little as 
10 from the end-on view. From the clear indications 
of radial motions described below, we can be certain 
that the bar is at least partially end-on, although it is 
difficult to gauge how much. We explore this in our 
models. 
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Fig. 6b. 



First we examine a line of sight or end-on orientation 
and assume that the half-length of the bar roughly ex- 
tends to 1.25' (4 kpc), or just within a ring of higher- 
density in the disk (Figure [5]). A linearly increasing ro- 
tation curve is used, which matches the rest of the disk 
at radii larger than 45". (Note: aside from the bar, no 
H I is included in the main disk interior to this radius.) 
Adding a bar in this way allows for some spread in the 
velocity range (Figure ITOl column 3), while remaining in 
a narrow range along the minor axis. However, the full 
velocity spread seen in the data is not achieved. Such a 
model is improved by adding radial motions of ±20 km 
s _1 (Figure [101 column 4). Due to projection effects and 
limited spatial resolution, we cannot discern whether this 
is inflow or outflow. 

For ±30 ° and ±45 offsets from the line of sight (pos- 
itive corresponds to the near side of the bar being in the 
NW quadrant), we achieve a somewhat improved fit with 
radial motions of ±40 km s _1 . Finally, increasing the 
offset to ±60 ° from the line of sight requires increasing 
the radial motions to ±60 km s _1 . These results ignore 
any additional radial motions in the disk beyond the bar, 
which are shown to be necessary in § 14.1.51 These mo- 
tions in the disk render constraining the orientation and 
radial motions of the bar even more difficult. 

The length of the bar and the dispersion of the gaussian 



distortions remain constant in all cases. For the figures 
in this paper we display a —30° offset bar with a radial 
motion of —20 km s _1 . It should be noted that ±30° 
paired with +20 km s _1 produces the same results. The 
disk radial motions decrease the amplitude of the optimal 
bar radial motion, hence we choose —20 km s _1 for the 
latter instead of —40 km s _1 as mentioned above. 

4.1.4. Additional Morphological Modifications 

In the approaching half, a clump of brighter emission is 
seen between —20 and —50" above the plane of the disk, 
centered at a velocity of 1075 km s _1 (Figure El panel 
corresponding to 5'). This is not to be confused with 
the slant due to radial motions mentioned in § [31 and 
described in § 14.1.51 To sufficiently achieve both of these 
elements, an arc centered at —55 ° from the line of sight 
in azimuth within the disk, and spanning 45 , is added 
to the far quadrant of the approaching half in all models 
(Figure [9]). This is likely an azimuthal continuation of 
the larger radial extent of the receding half relative to 
the approaching half, or possibly spiral structure, rather 
than a unique feature. The modeling of such features 
as spiral structure will be explored in future HALOGAS 
papers. 

4.1.5. Radial Motions in the Disk 
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Fig. 7. — Channel maps showing the data, one-component (1C) model, and various flaring (F) models . T he fina l two columns (Fio and 
Fio + Lag) display the flare model with radial motions and then with the optimal lag described in ij |4.1.5l and l4.1.6l Only subtle differences 
are seen w ith the addition of a lag, all within 3-4' from the center, most notably in the bottom row. More convincing evidence is shown 
in Figure 1131 All models shown here include a bar with an orientation 30° from the line of sight and a radial motion of —20 km s™ 1 . 
Contours are as in Figure [2] Velocities are given in each panel of the first column. 
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Fig. 8. — Key parameters used in the optimal models, with the 
1C and F models differing only in scale height. These are rotational 
velocity (vrot), column density (Cdens), position angle (PA), incli- 
nation, and for the F model, scale height (SH). Inclination is shown 
using thickened lines. Note that the receding half extends further 
than the approaching half, which drops to zero in column density 
at approximately 30 kpc. Values for other parameters beyond this 
point are excluded in that half. 

There are several indications of radial motions beyond 
the bar. This is most readily seen in Figure [3] as well as 
in Figure [101 Note first the slope of the highest contours 
on the systemic side of the approaching half, particularly 
evident in the panels corresponding to 0' and 1'. This is 
appropriately reversed for the receding half. These are 
indicative of radial inflow at small to intermediate radii. 
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Fig. 9. — A face-on view of the optimal model as seen from 
a viewer located in the SW quadrant. Note the extent of each 
half, the bar at the center with a half-length of 1.25', the ring 
of higher density at approximately 1.7', as well as the additional 
arc added to the far side of the approaching half. The black dot 
indicates the original orientation of the observer. Contours begin 
at 7.7 X 10 15 cm -2 and increase by factors of 2. Grid units are in 
pixels and 1 pixel = 4". 

This is modeled by introducing an inflow starting with 
an extreme value of — 50±5 km s _1 at 1.75' (5.5 kpc) 
and quickly decreasing with radius before becoming un- 
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Fig. 10. — Minor axis position-velocity diagrams showing models with no bar, a bar oriented 30° from the line of sight, the same bar 
with an inflow of —40 km s _1 , and the same bar with an inflow of —20 km s — 1 as well as additional radial motions that are not associated 
with the bar. These additional radial motions are described in § 14.1.51 Note the large spread in velocities in each panel related to the bar 
itself, as well as the slanting of the contours at 1175-1250 km s — 1 , indicative of additional radial motions. Contours are as in Figure [2] 
Dashed lines indicate the major axis and systemic velocity. 

be seen in Figure QTJ 

It is difficult to tell with certainty if the observed ra- 
dial motions are due to streaming along spiral arms or 
another source. However, the substantial inflow located 
within and just outside the ring of higher density in- 
dicates that the inflow in this region may be true ax- 
isymmetric inflow or associated with the bar rather than 
spiral arms. 



detectable beyond 2.75' (8.6 kpc). This inflow is symmet- 
ric in both halves and is very likely associated with the 
ring of higher density. The model in the fourth column 
of Figure [3] (Fi) includes all of these motions. Arrows in 
Figure indicate the significant but often subtle effects 
of adding such inflow. The directions of the arrows corre- 
spond to the direction in which emission is displaced due 
to radial motions. Additional isolated regions of inflow 
are seen in each half, but these are not symmetric and 
are likely associated with spiral structure. It should be 
noted that a more modest inflow of 10-15 km s _1 spread 
over a larger radial range could not reproduce what is 
seen in the data. 

In addition to the indications of radial inflow in the 
inner parts of the galaxy, indications of outflow are seen 
at large radii. This is best seen as the overall negative 
slope on the systemic side, opposite in direction from that 
produced by the inflow, ±40-80" off the plane of the disk 
in panels between -5' to +5' (once again indicated by 
arrows). This is modeled by allowing a radial outflow of 
10 km s^ 1 beginning at 6.75' (21.2 kpc), and extending 
outward to the edge of the disk in each half [Figure G2 
column 5 (Fio)]. 

The distribution of radial motions in the models may 



4.1.6. The Addition of a Lag 



As previously discussed in § 14.1. H and seen in the bv di- 
agrams in Figure [3j the data are relatively thin near the 
terminal side compared to the systemic. This is some- 
what remedied by the addition of a flare (§ 14.1.21) . but 
improvements may still be made through adding changes 
in rotational velocity with height above the disk. This 
causes observed velocities of emission at high z to be 
displaced towards the systemic side in the bv diagrams, 
resulting in a thinning in the width near the terminal side 
and a thickening on the systemic side. This is best seen 
in the bv diagrams closer to the center, shown in Fig- 
ures I12al and b where the contours on the terminal side 
become more rounded to better match the data with the 
addition of a lag. This is also evident at the terminal 
edge of lv diagrams, where emission is displaced to lower 
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Fig. 11. — Radial motions in each half, excluding those associ- 
ated with the bar. Negative values indicate inflows. Due to limited 
resolution, and complications from the bar, no radial motions are 
constrained for radii less than 45" (2.4 kpc). The sharp dip coin- 
cides with the ring of higher density, while the shallower dips are 
likely associated with spiral structure. 

velocities. It should be noted that these effects arc subtle 
but cannot be reproduced by changing the flare param- 
eters. If a lag is present in the data, then we should see 
some or all of these features, provided there is sufficient 
resolution. However, given how thin and far away the 
disk is, our resolution is lacking, so our error bars are 
large. We find an optimal lag peaking at -40 ijo km 
s _1 kpc -1 between 1.25' and 4.75' (3.9 and 14.9 kpc) in 
the approaching half, and -30 km s _1 kpc -1 between 
1.25' and 4.25' (3.9 and 13.4 kpc) in the receding half. 
These quickly decrease in magnitude within the speci- 
fied range in both halves, as shown in Figure 1131 This 
lag distribution is again represented in Figure [TJ] as an 
azimuthal velocity that rises more steeply at high z. 

For completeness we also include a model with a con- 
stant global lag throughout the galaxy in the sixth col- 
umn of Figure [3] Such a lag overestimates the rounding 
on the terminal side in panels corresponding to ±5' and 
±7', emphasizing the need for radial variation in the lag. 

4.2. On the Uniqueness of the Models 

A rigorous statistical analysis of the fitting of tilted 
ring mod els is generally om it ted from the litera- 
ture (e.g. iBarbieri et all 120051 lOosterloo et~aT1 120071 
iRand fc Benjaminll2008t ). Such studies employ a simple 
model focusing on the quantification of certain features 
in observed galactic disks, driven by a specific scientific 
question. In our case the main aim is to constrain the 
vertical structure of the H I disk, which means giving sig- 
nificant weight to faint emission. If fitting were fully au- 
tomated, a goodness-of-fit criterion such as a % 2 (which 
TiRiFiC provides) could be used to distinguish mod- 
els by significantly downweighting the bright emission 
that would otherwise dominate the statistic and render 
it meaningless for our purposes. In principle, a reduced 
X 2 with such a weighting could also account for the vary- 
ing number of free parameters among models. However, 



it is generally much more efficient to use the judgment 
of the modeler and adjust certain parameters manually. 
But this makes an accurate accounting of the degrees of 
freedom in each model difficult. For instance, the ra- 
dial scale height dependence (Figure [8]), shows eight val- 
ues and eight inflections, indicating sixteen values were 
varied. However, for this and other parameters, finer 
variations may have been considered and discarded as 
providing insignificant improvement. Likewise, initially 
most parameters for the approaching and receding mod- 
els are the same, and then some parameters varied as 
needed, but some deviations between halves may have 
been tested and discarded to keep the models simple. If 
so, the number of parameters allowed to vary was larger 
than the final model would suggest. 

We therefore rely on fitting by eye to guide our choice 
of models. A visual comparison of data and models pro- 
vides a powerful method for deciding how models must 
be altered to fit the data, easily allowing attention to be 
focused on regions of interest. One statistic we could ex- 
amine, regardless of the number of degrees of freedom, 
is the rms of the residual cube, downweighting regions of 
bright emission. In fact, when we examine this statistic, 
using only emission above 3er, we indeed find a significant 
7% reduction for models which include a flare (all models 
beginning with "F" ) rather than a constant scale height 
(1C). However, among the flaring models, the numerical 
difference due to the addition of radial motions, a bar, or 
a lag is negligible. Ne verth eless, improvement is clearly 
seen in Figures |3l QUI Il2bl respectively. In spite of such 
limitations, we address the need for such features here. 

Firstly, the additional arc in the receding half f§ l4.1.4p . 
which is included in all models presented in this paper, 
can be justified rather simply. There are clear indications 
that a feature of higher flux density exists in this region, 
but only in one quadrant of the disk (Figure El panel 
corresponding to 5'). The feature could be reproduced 
with an overall change in position angle, allowing flux 
to be redistributed at high z, but as this feature is seen 
only in a select region of the disk, such a change would 
worsen the fit in other parts of the disk. Additionally, 
a change in position angle for only this region would be 
unphysical and was seen to be a less desirable match 
during the modeling process (not shown). 

The addition of radial motions duplicates slants, 
slopes, and curves in the data (Figure [3]) that are other- 
wise not reproducible. One could argue for changes in the 
position angle near the center to reproduce these slants, 
but through modeling, we see that such a change would 
not fit the data as well. Furthermore, warping near the 
center, as well as both the positive and negative changes 
in position angle required to reproduce the correct slant- 
ing would be unc haracteristic of warping typically seen 
in other galaxies ([Briggs 1993). 

It is clear that a bar exists in NGC 4565, 
both from previous works at other wavelen gths 
(jKormendv fc Barentine1l2010L iNeininger et al.lll996ft . as 
well as our models. Substantial improvement can be 
seen in Figure QUI although due to projection effects, the 
uniqueness of this bar is difficult to constrain. Thus we 
provide generous uncertainties for the quantities involved 

(§ra. 

Finally, the lag is best constrained using faint emis- 
sion, high above the plane of the disk. For this reason, 
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Figure [T2b] 
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Fig. 12b. — Magni fied panels corresponding to those with rect- 
angles in Figure I12al Note the rounding and slope of the contours 
near the terminal side in each panel. 
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Fig. 13. — The optimal, maximum and minimum lag distribu- 
tions in each half. Due to limited spatial resolution, no lag is 
constrained for radii less than 45" (2.4 kpc). 
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Fig. 14. — Azimuthal velocities in the midplane and 20" (~1 kpc) 
above the plane of the disk for the approaching (a) and receding 
(r) halves. 

improvements due to a lag do not exhibit substantial 
changes in the rms of the residuals. Nonetheless, sub- 
tle improvements are seen in Figure |6a] and b, as well 
as Figure [l"2"al and b. The characteristics of the lag may 
initially be interpreted as H I that appears (at least in 
projection) above the plane of the disk that is rotating 
more slowly than H I in the midplane. The issue of pro- 
jection effects could hinder the uniqueness of our models 
if such effects arc not considered carefully. However, sig- 
natures in both bv diagrams and channel maps described 
at the beginning of § 14.11 allow us to successfully rule out 
inclination effects in this case. From this we may con- 
clude that the observed effects are extremely likely due 



to a lag. A range of possible values for this lag is given 
in Figure [T3l 

5. THE COMPANION GALAXIES 

While NGC 4562 is not our primary target, some infor- 
mation may be extracted from a basic model. We assume 
a distance of 12.5 Mpc based on values obtained from the 
NASA/IPAC Extragalactic Database (NED). The total 
Hi mass for this galaxy (including a factor of 1.36 to 
correct for He) is 3.5x10 s M Q . Our models include flux 
extending out to a radius of 2.75' (10 kpc). This small 
angular size limits our modeling capabilities in this case. 
The scale height is found to be about 7" ( 400 pc). The 
galaxy is slightly inclined at 82 . There are no indica- 
tions of a strong warp component along the line of sight, 
although there exists a slight warp component across it. 
The rotation curve is relatively flat, and peaks at 65 km 

No model is created for IC 3571. However, assuming 
approximately the same distance as NGC 4565, its Hi 
mass is found to be 5.6 xlO 7 M & and the extent of its 
largest axis is approximately 1.8' (5.8 kpc). 

6. DISCUSSION 

6.1. The Warp in NGC 4565 and Interactions With 
Companions 

The warp in NGC 4565 is highly asymmetric, not only 
extending further and reaching higher values in position 
angle on the receding half, but also flattening at large 
radii on that side as well (Figure [5J . There was a hint 
of thi s flattening in a few channels as noted by IRupenl 
(|1991f ). most prominent in the channel corresponding to 
Vhd = 1459 km s _1 . Here we see it from Vhei = 1440 to 
1473 km s" 1 . A s imilar flattening is also seen in NGC 
5907 (|Sancisilll983l) . 

Additionally, there is a slight rise in inclination for radii 
larger than 8.75' (27.5 kpc), causing the galaxy to be- 
come more edge-on. This is only seen in radii beyond 
the extent of the emission in the approaching half, thus 
it is only seen in the receding half. 

It is difficult to determine the exact cause of the asym- 
metries in the warp, but it is possible that the companion 
galaxies play some role. NGC 4565 is clearly interacting 
with IC 3571, as can be seen in Figured] by the bridge of 
material between them. Additionally, material closer to 
the plane of the disk appears to be displaced above the 
plane, toward IC 3571 (Figure [T]). Given the location 
of IC 3571, it may contribute to the asymmetry in the 
component of the warp across the line of sigh t. This was 
also noted bv Ivan der Hulst fc Sancisil ([20051 ). 

Emission likely tied to IC 3571 as part of the bridge be- 
tween the two galaxies could not be reproduced in mod- 
els, which is unsurprising. The velocity range associated 
with this emission is reasonably close to that of the main 
disk of NGC 4565 (Figure [15]), and is likely only a dis- 
turbance of gas within the disk as opposed to accretion 
from IC 3571. 

6.2. The Flaring Gas Layer 

The exponential scale height in NGC 4565 begins to 
increase in both halves from an initial 4" ( 200 pc) near 
a radius of 3 arcminutes (9.4 kpc), eventually reaching a 
value of 12" ( 600 pc) at a radius of 7 arcminutes (22.8 
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Fig. 15. — Position-velocity diagrams along the minor axis show- 
ing the bridge between NGC 4565 and IC 3571 in the smoothed 
(44" X 34") cube. Contours begin at 2a (0.48 mjy bm" 1 ) and 
increase by factors of 2. A — 2cr contour is shown as dotted lines. 
Major axis offsets are given at the top of each panel. 

kpc). This is just within the optical radius (8.1' or 25.5 
kpc). 

Fl ares are seen in othe r galaxies, including M 
31 dBrinks fc Burton! 1 19841) and the Milky Way 
(jKalberla et al.l 120071 and references therein). However, 
there exist galaxies for which extensive H I modeling has 
been done that do not show flaring ga s layers, includ- 
ing NGC 5746 pand fc Beniaminl l2008l) and NGC 4559 
(jBarbieri et al.ll2005[ ). indicating that flares are common 
although not ubiquitous. 

6.3. Halo Trends 

Our models yield no evidence for a massive, extended 
H I halo in NGC 4565. This is shown by the vertical pro- 
file in Figure [5] where there are no indications of a second, 
extended component that exists globally throughout the 
data. Furthermore, a second component is not seen to 
improve position-velocity diagrams or channel maps (not 
shown). 

We return to the issue discussed in §[T]of whether there 
is any connection between H I halos and disk star forma- 
tion activity. It is difficult to extract any concrete trends 
regarding H I halos with such a small number of galax- 
ies having been observed to this depth and subsequently 
modeled in detail. The situation is in the process of be- 
ing remedied via the HALOGAS survey and Expanded 
Very Large Array (EVLA) observations involving some 
HALOGAS team members. In Table[3]we compare prop- 
erties including total H I mass, total infrared luminosity 
per unit area, scale heights, and lags of NGC 4565 with 
those of previously modeled edge-on galaxies. 

The total infrared luminosity per unit area for each 
galaxy is calculated in this work unless otherwise stated. 
Th ese calculat i ons ar e done based on equations found 
in iDale et al.l (I2009D. and involve the optical area 
(|de Vaucouleurs et al.lll99ll in most cases; additional ref- 
erences are provided in the table). Spitzer MIPS data 
are used if available for a given galaxy, otherwise IRAS 
data are used. 

Immediately apparent is the pattern of higher infrared 
luminosity per unit area (i.e. star formation rate per unit 
area) and the presence of extra-planar H I in the form of 
a halo. The obvious exception is UGC 7321. NGC 5775 
could be considered an exception, since a global H I halo 
is absent. However, in that case, prominent extra-planar 
features including extensions and arcs are present. 

Still, we must emphasize that this result is based on 
only a few galaxies and modeling methods differ amongst 
them. Further analysis of the HALOGAS sample, in 



which individual galaxies are modeled in a fashion con- 
sistent from galaxy to galaxy, will provide a stronger test 
of this trend. 

6.4. The Lag in NGC 4565 and Other Galaxies 

We do not detect substantial extra-planar H I in the 
form of a global halo, which indicates that there is no 
strong disk-halo cycling in NGC 4565. However, the Hi 
has some vertical velocity dispersion and experiences a 
lag when it rises above the plane of the disk. The lag we 
detect is within the innermost 4.75' (14.9 kpc) and 4.25' 
(13.4 kpc) in the approaching and rec eding halves, re- 
spect ively. As is the case for UGC 1281 (iKamphuis et al.l 
I20T1 and NGC 4244 (jZschaechner et al.l 1201 1[ ). this in- 
dicates that substantial extra-planar gas need not be 
present in order for a lag to be observed. (Although 
for UGC 1281, a larger line of sight warp component is 
favored over a lag, yet a lag cannot be ruled out.) 

The H I lag detected in NG C 4565 is ver y steep . This 
would be consistent with the iHeald et al.l (|2007| ) trend 
in which DIG lags are seen to steepen with decreasing 
SFR, provided that H I and DIG lags share a common 
origin. Although the DIG lag in NGC 4565 is unknown, 
modeling of the DIG within the HALOGAS project is 
currently underway (Wu et al. 2013, in prep). 

A radially varying lag is seen in NGC 4565, which we 
will now discuss. The onset of this shallowing is at 2.75' 
(8.6 kpc) in the approaching half and 2.25' (7 kpc) in the 
receding. Given the large uncertainty associated with the 
lag in general, it is difficult to constrain the manner in 
which it shallows, although outside of 4.75' (14.9 kpc) in 
the approaching half, and 4.25' (13.4 kpc) in the receding 
half, it appears to be zero. 

An entirely geometric interpretation of the shallowing 
of a lag at large radii, without any consideration of ef- 
fects such as extra-planar pressure gradients, drag, or 
magnetic fields, considers the relationship between the 
rotational velocity and the gravitational potential: gas 
rotating high above the plane of the disk is further from 
the gravitational center than the gas at the same ra- 
dius in the midplane, resulting in a slower rotational ve- 
locity. This effect is more pronounced at smaller radii, 
where a certain height above the disk is comparable to 
the radial distance from the center, resulting in a larger 
relative change in the gravitational potential, and thus 
a steeper lag . This effect may be seen in Figure 3 of 
iCollins et al.1 (|2002f ) and may in part be why a lag is 
only seen in the inner third of the disk in NGC 4565. 
However, it should again be noted that such purely ballis- 
tic models have under-predicted observed lag magnitudes 
(|Fraternali fc Binnevll2006llHeald et~aT1l2007l) . indicating 
that additional factors must be considered. 

This radial variation could also be affected by an extra- 
planar pressure gradient if the cycling gas does not orbit 
purely ballistically a nd is subje c t to h ydrodynamic influ- 
ences. According to iBeniaminl (|2002l ). a pressure gradi- 
ent directed radially inward may cause a steeper global 
lag. An outward pressure gradient may result in a shal- 
lower global lag, or even an increase in rotation velocity. 
If such pressure gradients vary radially themselves, then 
a radial change in the lag may occur. Given how little 
is currently known of halo pressure gradients, as well as 
the potentially substantial impacts even small pressure 
gradients may have, it is difficult to say to what degree 
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TABLE 3 

Disk and Halo Properties of Edge-on Galaxies 



Gcilcixy 


H I mass 




1 p XJ a 


Thickest Component 


H I Lag 


rxei. 




10 9 M Q 


10 40 ergs -1 kpc 2 


or F 


Scale Height (kpc) 


km s~ 1 kpc~ 1 




UGC 1281 


2.75 b 


n/a c 


F 


0.18-^0.31 


n/a 


1, 2 


UGC 7321 


1.1 


0.35 


H, F 


2.2 d 


< | - 25| 


3 - 7 


NGC 4244 


2.5 


0.53 


1C 


0.56 


-9 


3, 5, 8, 9 


NGC 5746 


10 


1.16 


IC 


0.4 


n/a 


3, 4, 10, 11 


NGC 4565 


9.9 


1.45 


F 


0.2 ->0.6 


-30, -40 


3, 9, 11, 12 


Milky Way 


8 


3.0 


H, F 


1.6 


-22,-15 


11, 13, 14 


NGC 3044 


3.0 


6.62 


H 


0.42 d 


n/a 


3, 4, 14, 15 


NGC 891 


4.1 


8.71 


H, F 


1.25->2.5 


-15 


3, 8, 12, 16 


NGC 5775 


9.1 


25.6 


IC + features 


1 


n/a 


3, 4, 14, 17 



a 1-component (IC), halo (H), flare (F), or a combination of these. 

(Private communication) P. Kamphuis. 
c Appropriate Ltir is not available for this galaxy. However, based on Ha, its star formation rate is 0.0084 Mq yr — 1 
( van Zee 20011), rendering it among the lowest of the galaxies presented here based on calculated rates for the rest. 
d Gaussian as opposed to exponential layer. In the case of NGC 891, sec Oostcrloo ct al. 2007 for details on the distribution 
used to fit the halo component. 

References: (1) Kamphuis et al. 2011, (2) van Zee 2001 , (3) This work, (4) do Vaucouleurs ct al. 1991, (5) Dale et al. 
2009, (6) Uson & Matthews 2003, (7) Matthews & Wood 2003, (8) Heald et al. 2012, (9) Zschaechner ct al. 2011, (10) 
Moshir et al. 1990, (11) Rand & Benjamin 2008, (12) Rice et al. 1988, (13) Marasco & Fraternali 2011, (14) Sanders et 
al. 2003, (15) Lee & Irwin 1997, (16) Oostcrloo et al. 2007, (17) Irwin 1994 



these will alter the characteristics of any given lag. 

In Hi, radia lly shallowing la g s ar e also seen 
in N GC 4244 dZschaechner et al.l [201lh and NGC 
891 (IQosterloo et al l 120071) . and the Milky Way 
iMarasco fc Fraternali! (|2011f l. In NGC 4244 the lag de- 
creases from —9 km s _1 kpc -1 to —5 km s _1 kpc -1 and 
—4 km s _1 kpc -1 in the approaching and receding halves 
respectively. This would correspond to a shallowing of 
2-2.5 km s~ x kpc~ 2 (although, as stated in that paper, 
and is the case for NGC 4565, the uncertainties in the 
radial variation are large). In NGC 891, the gradient in 
the inner regions is —43 km s _1 kpc -1 and decreases to 
— 14 km s~ x kpc -1 in outer radii, corresponding to a shal- 
lowing of 2.5 km s _1 kpc -2 . In the Milky way, the lag 
starts with a magnitude of approximately —45 km s -1 
kpc -1 and reaches a value of —15 km s _1 kpc -1 near 
4.5 kpc, resulting in a more rapid shallowing of 6.7 km 
s" 1 kpc -2 . The lag in NGC 4565 also decreases rapidly 
with radius. Considering the onset of the shallowing of 
the lag in each half and their respective lag magnitudes 
(Figure [T3f . this corresponds to a shallowing of 6.4 km 
s" 1 kpc~ 2 and 4.8 km s _1 kpc -2 in each half for our op- 
timal models. Judging by the minimum and maximum 
lag models, the shallowing could be between 6.4 and 38 
km s~ x kpc -2 in the approaching half (although the lat- 
ter value is very unlikely as it would imply that a very 
rapid drop of —60 km s -1 kpc -1 occurs in a very nar- 
row annulus of 1.6 kpc, which does not seem physically 
realistic) and between 4 and 7.6 km s _1 kpc -2 in the 
receding. Thus, the lag in NGC 4565 appears to shallow 
more rapidly than in NGC 4244 and NGC 891, but at a 
similar rate to that of the Milky Way. 

In NGC 4244, the onset of the shallowing of the lag 
is seen just within the optical radius, whereas the shal- 
lowing of the lag in NGC 891 is closer to its center. In 
the Milky Way, the lag begins to shallow immediately at 
the center, and reaches a constant value near 4.5 kpc. In 
NGC 4565, the shallowing of the lag begins at approxi- 
mately 2.5' (7.9 kpc), which is less than half the optical 
radius of 8.1' (25.5 kpc). The lag then reaches a con- 



stant value of zero near 5' (15.7 kpc). Thus, there does 
not yet appear to be a characteristic radius at which the 
lags begin to shallow, or where the lags reach constant 
values. 

In addition to these edge-on examples, NGC 4559, a 
moderately inclined ga laxy also shows evide nce for a ra- 
dially decreasing lag (|Barbieri et al.l 120051) . Although 
a true velocity gradient in height abov e the plane of 
the di sk cannot be measured in this case, iBarbieri et al.l 
(2005) find a decrease in the rotational velocity of the 
thick halo component compared to the thin disk of 
roughly —60 km s _1 near the center of the galaxy, and 
—20 km s _1 at large radii. As long as the vertical ex- 
tent of the lagging halo does not fall in a similar way, 
this drop would indicate a true radially decreasing lag. 
However, because of its inclination, NGC 4559 does not 
yield an unambiguous case of a radially varying lag. 

This small sample shows evidence for discrepancies in 
the nature of the radial variations in lags in each galaxy, 
including differences in the radius of the onset of shallow- 
ing and the rate of shallowing. Although three edge-on 
external galaxies as well as the Milky Way clearly show 
a shallowing lag with increasing radius, a larger sample 
must be studied in order to extract reliable trends. 

6.5. The Bar in NGC 4565 

We hnd kinematic evidence for radial motions, which 
we associate with a bar. Given the inclination of NGC 
4565, we only observe a rather limited 2-D projection; 
not nearly enough to determine its length and thickness 
using only spatial information, so we must rely on its 
kinematics. However, observations at other wavelengths 
are considered in order to help constrain morphological 
properties. In the optical, the bar is seen as a boxy bulge 
(jKormendv fe B arcntinc 2010). Additionally, a ring may 
be seen in higher resolution 24/Lt MIPS (NASA/IPAC In- 
frared Science Archive) data (Fig ure [T6"|). as well as CO 
and 1.2 mm continuum emission ()Neininger et all 119961 
Yim et al. 2012 in prep). The radius of the ring (1' to 
1.7' in the 24fi data, and peaking from 1' to 1.5' with an 
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Fig. 16. — Grayscale of 24 MIPS data with Hi contours 
overlaid starting at 3.0 X 10 19 cm~ 2 and increasing by factors of 2. 
Note the prominent ring with a radius of just under 2 arcminutes. 

outer radius of 3' in the CO data) is consistent with the 
ring of higher density seen in H I, and is thus consistent 
with the length of the bar we have modeled. 

Parameters related to the bar are difficult to constrain, 
thus we will only state what may be reliably extracted 
from these data. We can say that the bar exists with a 
half length of approxim ately 1.25' (4 kpc). According to 
iNeininger et al.l (|1996f) the bar's orientation is between 
and ±45 from end-on. Due to limited resolution 
and projection effects, we cannot improve upon this. In 
fact, by adjusting the rotation curve along the bar and 
radial motions within the disk, these data allow for an 
even larger azimuthal range. 

There appear to be radial motions up to ±60 km s _1 
within the bar. For these radial motions we cannot dis- 
cern between inflow and outflow, largely because we can- 
not determine the sign of the orientation due to projec- 
tion effects, although wi thin the length of th e bar, a net 
inflow is expected (e.g. lAthanassoulal [19 92') . Further- 
more, even without changing the sign of the orientation, 
changing the direction of radial motions produces too lit- 
tle of an effect at most orientations to determine which 
is correct. This may be further complicated by radial 
motions that vary in magnitude along the bar, of which 
there is some indication in that further from the center 
there appears to be less need for such motions, but this 
is complicated by projection effects and is not enough to 
make meaningful constraints. These properties are ex- 
tremely interconnected, and our spatial resolution lim- 
ited, hence the large ranges presented here. Fortunately, 
as may be seen in Figure 1101 the effects from radial mo- 
tions outside the bar far outweigh any contributions from 
the bar itself. Thus, we must stress that the conclusions 
we draw from our final models that are not directly about 
the bar are unaffected by the bar itself. 

6.6. Radial Motions in the Disk 

Observations of purely axisymmetric inflows, not as- 
sociated with motions along a bar or sp iral structure, 
have been ambiguous (|Wong et al.l 12004] and references 
therein). Our work does not break this ambiguity, but 
rather provides additional constraints. 



We see clear indications of a net inflow beginning 
around the ring of higher density (1.75' or 5.5 kpc), and 
extending out to 2.75' (8.6 kpc). The extreme value of 
this inflow is — 50±5 km s _1 . Unlike the issues related 
to constraining the bar, the signatures for this inflow are 
clear and appear to be unique. There are additional re- 
gions in which the models are improved by inflow, but 
these are almost certainly due to spiral arms. Since this 
inflow is associated with the ring, it could be associated 
with the bar potential. Finally, there exists evidence for 
a net outflow of 10 km s _1 beginning near a radius of 
6.75' (21.2 kpc) and extending outward in both halves. 

The observed effects of the outflow are primarily seen 
above the plane of the disk at a height of 40 to 80" (2-4 
kpc). A range of inclinations corresponding to a warp 
component along the line of sight were tested, but none 
could replicate the asymmetric slant in the data. Fur- 
thermore, if the apparent offset from the plane of the disk 
were due to projection effects alone, then these would 
correspond to radii between 50 and 70 kpc, placing them 
well outside the disk. Therefore, these motions likely do 
exist above the plane of the disk, rather than being due 
to such effects. It is also possible that any outflow could 
be caused by material being pulled outward due to tidal 
interactions. 

Both NGC 4559 (IBarbieri et all 12001 an d NGC 2403 
(jSchaap et all 120001 iFraternali et al.l 120021 ) have shown 
similar indications of radial inflows, estimated between 
10 and 20 km s _1 . In each case, the authors note that the 
inflow is high above the plane of the disk. As mentioned 
above, our data and models indicate that most of the 
outflow in NGC 4565 is in the thicker, flaring layer. The 
inflow however, is close to the plane of the disk and is 
significantly higher than these values. 

7. SUMMARY AND CONCLUSIONS 

Through our models we determine the following: 

1. In NGC 4565 we see a flaring layer, but no massive, 
extended extra-planar H I layer throughout the galaxy. 
There are two companion galaxies, one of which is clearly 
interacting with NGC 4565 via a bridge of extra-planar 
material. 

2. We detect a lag peaking in magnitude at -40 t 2 o 
km s _1 between 1.25 and 4.75' (3.9 and 14.9 kpc) on the 
approaching half, and -30 km s~ x between 1.25 and 
4.25' (3.9 and 13.4 kpc) on the receding half. Both lags 
quickly decrease with radius and are zero at radii larger 
than the upper end of the specified range. 

3. A bar was previously known to exist within NGC 
4565 for which we model a half-length of 1.25' (4 kpc). 
Our models indicate radial motions along this bar of 20 
to 45 km s _1 . A ring of relatively high density is seen 
near the end of the bar. 

4. Our models indicate a radial inflow beyond the bar 
that may be associated with the ring of higher density. 
This inflow starts at a radius of 1.75' (5.5 kpc) with a 
value of — 50±5 km s _1 and decreases with radius before 
becoming undetectable outside of 2.75' (8.6 kpc). A net 
outflow of 10 km s _1 may be present at large radii, most 
prominently seen 2-3 kpc above the plane of the disk. 
The origin of this outflow is unclear. 

At this time, we are still limited to a small number of 
galaxies for which deep observations exist and detailed 
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models have been made. Even fewer of these, such as 
NGC 4565, are edge-on. However, in edge-ons, a trend is 
now seen suggesting some connection between the pres- 
ence of substantial extra-planar H I and star formation 
rate per unit area within a given galaxy. Additionally, we 
now see radially shallowing lags in three edge-on galax- 
ies, and possibly one that is moderately inclined. Further 
work with the full HALOGAS sample will either dismiss 
or strengthen these early results. 
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